Magnetodielectric and magnetoelastic coupling in TbFes (603)4 
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We have studied the magnetodielectric and magnetoelastic coupling in TbFes (603)4 single crys- 
tals by means of capacitance, magnetostriction and Raman spectroscopy measurements. The data 
reveal strong magnetic field effects on the dielectric constant and on the macroscopic sample length 
which are associated to long range magnetic ordering and a field-driven metamagnetic transition. 
We discuss the coupling of the dielectric, structural, and magnetic order parameters and attribute 
the origin of the magnetodielectric coupling to phonon mode shifts according to the Lyddane- Sachs- 
Teller (LST) relation. 

PACS numbers: 77.80.-c, 61.10.Nz, 77.84.-s 
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I. INTRODUCTION 

Many novel complex systems are investigated in 
the research on magnetoelastic and magnetodielectric 
couplingir— . Among them, rare-earth iron borates 
RFe3(B03)4 constitute a family of non-centrosymmetrric 
oxides in which 4f-moments are embedded in a network 
of id Fe spins. These compounds crystallize in trigonal 
R32 space group at high temperatures^!^. The members 
of the family having rare-earth elements with ionic radii 
smaller than Sm undergo a structural phase transition 
lowering the crystal symmetry from R32 to P3i21<^ii. 
Based on single crystal diffraction studies, it has been 
proposed that the structural transition results in an an- 
tiferroelectric stated. Magnetic exchange interactions for 
R = (magnetic) Rare Earths are complicated due to the 
4f-3d interactions and the anisotropy of the 4f-moments 
which, e.g., affects the direction of Fe-spin in the ordered 
phase. 

Structurally, TbFe3(B03)4 exhibits edge sharing FeOg 
polyhedra which are forming helicodial chains along the 
c axisi. Fe spins order antiferromagnetically at Tn ~ 39 
K along the helicoidal chains of Fe atoms with the spins 
parallel to the hexagonal c axis while Tb moments are 
polarized from Fe ordering^"—. When a magnetic field 
is applied along the easy axis c, a field induced mag- 
netic transition is observed at temperatures below Tn, 
i.e. there is a spin-flop of the Fe spins while the Tb mo- 
ments become fully aligned. 

So far R = Gd, Pr, Nd and Tb compounds have 
been studied in the context of multiferroics and mag- 
netoelectrics. In fact, the space groups of both high and 
low temperature phases are not polar despite being non- 
centrosymmetric, therefore the family cannot be consid- 
ered as multiferroic. On the other hand, it has been 



shown for R — Gd, Pr, Nd and Tb that magnetic field 
induces electrical polarization, marking the materials as 
magnetoelectrioii"— . Also, sizeable spin-lattice coupling 
was demonstrated for some of these compounds using 
magnetostriction measurementsi^^. Within the family, 
the dielectric properties were reported only for R = Gd 
where the dielectric constant was found to decrease with 
decreasing temperature. The onset of magnetic ordering 
causes an anisotropic anomaly in e. In the magnetically 
ordered phase Gd 4f - Fe 3d magnetic interactions cause 
a spin reorientation transition which is coupled to the 
dielectric constant via spin-lattice interaction»iii^. Ap- 
plication of magnetic field strongly shifts the spin reori- 
entation transition, giving rise to a 1% magnetodielectric 
effeclfi^. Here, we present a detailed experimental study 
on the magnetodielectric and magnetoelastic coupling in 
TbFc3(B03)4. We have measured the dielectric constant 
and the magnetostriction and we show by analysis of Ra- 
man spectroscopy data that the magnetodielectric cou- 
pling is mediated by strain and caused by shifts in the 
corresponding transverse optical (TO) phonon modes via 
Lyddane-Sachs- Teller (LST) relationship. 



II. EXPERIMENTAL 

The single crystals we used are from the same batch 
as in Ref. |^. The capacitance of the sample was mea- 
sured using a home-made insert and Andeen-Hagerling 
AII2500A bridge in a 18 Tesla Oxford cryostat between 
5 and 300 K. The measurement frequency was 1 kHz. 
Magnetostriction measurements were carried out with a 
capacitance dilatometer— in the same cryostat and with 
the same capacitance bridge. Our dilatometer applies 
the tilted plate principle with the sample placed in the 



open center of a ring-like capacitor made from silver. It 
enables to follow the length changes of the sample at 
constant temperature upon changing the external mag- 
netic field. The sample length was 1.23 mm. The field 
was changed quasi-statically from up to 14 T and back 
with a rate of 0.01 T/s. The length changes of the sam- 
ple arc calculated from the capacitance changes measured 
by a temperature stabilized capacitance bridge Andeen- 
Hagerling with a resolution of 5*10 ~^ pF. Thus, length 
changes of less than 0.01 A can principally be resolved. 
Due to mechanical vibrations, etc., the resolution is lim- 
ited to 0.1-1 A in practice. The measurement frequency 
was 1 kHz. The field dependence of the magnetization up 
to 15 T was studied in a home-built VSM magnetome- 
ter ,i^ 

The Raman measurements were performed in a 
backscattering configuration, using a three-grating 
micro-Raman spectrometer (T64000 Jobin Yvon) 
equipped with a liquid nitrogen cooled charged coupled 
device (CCD) detector. The frequency resolution was 
better than 2 cm~^ for the frequency region considered. 
The triple grating configuration allows the analysis of a 
broad spectrum from 10000 cm~^ down to 4 cm~^. The 
sample was placed in an optical microscope cryostat. 
The temperature was varied from 4 to 55 K, with a sta- 
bility of 0.1 K. The scattering was excited by the second 
harmonic light of a Nd:YV04 laser (532 nm), focused 
to 50 /im^ with the power density on the sample kept 
below 0.1 mW//im^. The polarization was controlled 
both on the incoming and outgoing beams giving access 
to all relevant polarization combinations. The Raman 
measurements reported in the paper are performed with 
both the light's k- vector and the polarization of the 
incoming and outgoing field in the BO3 planes. 
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FIG. 1: Temperature dependence of the dielectric constant 
of TbFe3(B03)4, for E \\ c (a) and £; _L c (b). The insets 
highlight the data at low temperatures. The red line show 
fits according to Barrett's function (see the text). 



III. RESULTS AND DISCUSSION 

Fig. 1 (a) shows the temperature dependence of the di- 
electric constant measured along the c axis. In the whole 
temperature regime under study, ec decreases upon cool- 
ing. In addition, we observe an abrupt decrease at the 
structural phase transition at Ts = 202 K. Also the mag- 
netic ordering of the Fe spins at Tn « 39 K causes an 
anomaly in the dielectric constant (see the inset of Fig. 
1(a)) which is however much smaller than that at Tg. 
The temperature dependence of the dielectric constant 
measured perpendicular to the c axis in general resem- 
bles the behavior for E\\c but there are subtle differences 
at low temperature (Fig.[T]). We again observe a discon- 
tinuity at Ts but unlike the case i? || c there is only a 
very weak anomaly at Tn (cf. inset of Fig.[TJb)). In- 
stead, the data display a minimum above the magnetic 
ordering temperature, i.e. around 50 K, below which the 
dielectric constant increases. Upon further cooling, e^b 
decreases again below a broad maximum around 20 K. 

The effect of an external magnetic field on Ec applied 
parallel to the c-axis is shown in Fig.I^Ja). At T = 5 



K, the capacitance is rather field independent for small 
B\\c but a sharp increase appears around 4 T. At higher 
fields, the capacitance increases nearly linearly with a 
rather constant slope. The critical field of the jump in 
Ec increases upon heating and the size of the jump be- 
comes smaller. Above Tn, i-C in absence of long range 
magnetic order, the jump Ae vanishes. Note, that for all 
temperatures studied in the magnetically ordered phase 
the behavior in the high-field phase is nearly independent 
of temperature and the slope at high fields is very sim- 
ilar above Tn, too. The sharp anomaly in the dielectric 
constant is associated to a previously found field induced 
magnetic transition (FIMT) which is clearly visible in 
magnetization measurements (see. Fig.I^Jb)) 4 The mag- 
netization data show a huge jump of roughly 9 /le/f-u. at 
the FIMT, and similar to the jump in e the anomaly AM 
decreases and shifts to higher fields upon heating. The 
field induced inagnetic transition is not only associated 
to jumps in Cc and in the magnetization but it is also 
accompanied by strong magnetostrictive effects. As seen 
in in Fig. 2(c), the c-axis exhibits a jump-like decrease at 
the critical field. Again, the jump reduces upon heating 
and it is restricted to the magnetically ordered phase. In 



general, the c axis contracts in the whole field range stud- 
ied. In Fig. 2(d), the magnetic field dependence of eab is 
shown, with the magnetic field again applied along the 
c axis. We again observe that the magnetic field depen- 
dence is small in both the low-field and the high-field an- 
tiferromagnetic phase while the largest changes occur at 
the FIMT. However, unlike the measurement for Ec, the 
dielectric constant decreases with magnetic field in the 
field induced phase. The phase diagram shown in the in- 
set summarizes the transition fields of the FIMT vs. tem- 
perature as derived from dielectric constant, magnetiza- 
tion, magnetostriction and specific heat4. We note that 
sharp changes in magnetostriction and electrical polariza- 
tion at r = 4.2 K were reported recently in pulsed field 
studies, too.— i^ The reported anomalies in the magne- 
tostriction are however about 2 times smaller than found 
with our quasi-static set-up which might be associated 
to the different experimental approach, i.e. pulsed mag- 
netic field and usage of a piezoelectric sensor glued to 
the sample. In addition, there are also qualitative dif- 
ferences to our data, i.e. in ALa/La at B\\c > Be and 
pronounced hysteresis effects in the magnetostriction and 
the electrical polarization. 

To summarize the experimental data in Fig. [51 at the 
FIMT TbFea (603)4 displays a magnetically driven phase 
transition which is associated to sharp jumps in (1) the 
dielectric constant, (2) the magnetization, and (3) in the 
length of the c-axis. It has been shown earlier for the 
magnetic degrees of freedom that, at the FIMT, there is 
a spin-flop of the Fe spins accompanied by a magnetic 
flip of the Tb moments.— 

The pronounced magnetoelastic coupling and the re- 
sulting jump of the sample length at the FIMT might 
suggest that the observed step in the dielectric constant 
is only associated to the magnetostriction, i.e. 
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where eq is the permittivity of the free space, e is the 
dielectric constant, A is the area of the contacts, and d 
is the thickness of the sample. In general, a magnetic 
field induced jump in the sample length implies a capac- 
itance anomaly which can not be discussed in terms of 
the magnetodielectric effect. This scenario can however 
be clearly ruled out if the anisotropy of the magnetostric- 
tion and magnetocapacitance is considered (cf. Fig. [3]). 
If a magnetic field is applied along the c-axis, the length 
of the c-axis and the a-axis sharply shrinks and increases, 
respectively. Such an opposite sign of the length jumps 
is expected due to at least partial volume conservation of 
the magnetoelastic distortion. In contrast, both jumps 
in ec and e^b exhibit the same sign, i.e. magnetostriction 
has only a minor effect on the observed field dependence 
of e and cannot explain the changes in e{B). Therefore, 
the anomalies clearly point to a direct magnetodielectric 
coupling. 

As discussed above, the general behavior of the di- 
electric constant vs. temperature is rather isotropic 




FIG. 2: (a) The magnetic field dependence of the dielectric 
constant for E || c, B || c. (b) Magnetization (from Ref. [3), (c) 
Magnetostriction of the c axis, and (d) capacitance for E _L c, 
B II c. The inset displays tlie magnetic phase diagram derived 
from dielectric constant, magnetization, magnetostriction and 
thermal expansion studies. 
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FIG. 3: (a) Magnetocapacitance parallel and perpendicular 
to the c-axis, and (b) magnetostriction along the a- and c- 
axis, respectively, of TbFe3(B03)4, at T = 5 K, for magnetic 
fields applied along the c-axis. 



while there is a weak difference between ec(T) and ea(T) 
at low temperatures. Interestingly, despite the fact 
that in GdFe3(B03)4 the magnetic structure in the or- 
dered phase is different compared to TbFe3(B03)4, the 
anomaly at Tn as well as the low temperature behavior 
ic{T) and £ab{T) is roughly similar in both compounds.— 
In particular, both materials exhibit an increase of eab 
below Tn which in GdFe3(B03)4 is truncated by a spin- 
reorientation transition. For GdFe3 (303)4, the up- 
turn was suggested to originate from the polarization 
of Gd spins and the vicinity of a spin-reorientation 
transitioni^. In contrast, there is no spin-reorientation 
in TbFc3 (303)4 but as well the specific heat as the mag- 
netization data exhibit a broad Schottky-like anomaly at 
the temperature where there is the peak in Cab- Recently, 
the anomalies of the specific heat and the magnetization 
were explained in terms of temperature-driven popula- 
tion of the ground state of Tb ion split by the internal 
field of Fc spins— which suggests a similar, i.e. Schottky- 
like, scenario for the broad peak in in tab- 

In the following we will discuss the origin of the ob- 
served magnetodielectric coupling. In general, the di- 
electric constant is related to optical phonon frequencies 
via the Lyddane-Sachs- Teller relationship: 
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In this equation oj^ and lut a-re the long wavelength longi- 
tudinal and transverse optical-phonon mode frequencies, 
respectively, eg is the static dielectric constant, i.e. the 



dielectric constant at zero frequency, and £00 is the opti- 
cal dielectric constant. Following the early literature on 
BaMnF^, MnO^S and MnF2^, we have tried to relate 
the dielectric constant we measured to the relevant TO 
phonon modes. In order to do so, we have used a modi- 
fied Barrett equation as described in the aforementioned 
references as follows: 



e(r) = e(0) + A/[exp{hujo/kBT) - 1] 



(3) 



In this equation, ^ is a coupling constant and wq is the 
mean frequency of the final states in the lowest lying 
optical phonon branch. By fitting our data as indicated 
by the red lines in Fig.[T] we obtained A = 33.6, ujq ~ 
206 cm~i for £ab and A = 12.2, wq = 177 cm~^ for ec. 
We note the tranverse optical mode that we obtained 
from the Barrett function fits corresponds to an average 
of all the contributing transverse optical modes in the 
respective directional. 

We have done Raman spectroscopy measurements to 
look for direct evidence of spin-lattice coupling. Accord- 
ing to the previous report o^i^^ there are 59 transverse 
optical modes that correspond to eab i-c. phonon modes 
corresponding to a TO mode propagating perpendicular 
to the c axis. We have checked the temperature depen- 
dence of these modes and observed that only two low 
lying ones shift their frequency at low temperatures: the 
first one at 200 cm~^, the second one at 260 cm~^. The 
modes at other frequencies do not significantly change 
their frequency. Wc show in Fig.|4l[a) two representative 
low- lying phonon modes around 200 cm~^. Upon cool- 
ing below the magnetic ordering temperature, one of the 
modes with the frequency 206 cm~^ does not shift with 
changing temperature while the frequency of the lower 
lying transverse phonon mode at 200 cm~^ is shifting re- 
markably at Tn- Its shift provides direct evidence for the 
spin-lattice coupling since exchange striction stretches or 
elongates the bonds and hence shifts the phonon frequen- 
cies when magnetic ordering takes place. 

The anomalous softening of the phonon mode reported 
in Fig.HI^a) is a possible reason for the behavior of the di- 
electric constant. Below the antiferromagnetic ordering 
temperature Tn, the behavior of the dielectric constant 
perpendicular to c (as discussed in the previous section) 
is rationalized in the light of the Raman response. As 
suggested by the Lyddane-Sachs- Teller (LST) equation 
the dielectric constant is inversely proportional to the 
frequency of the TO-phonon mode, associated to it. We 
note that usually LO-modes and Eoo are assumed as tem- 
perature independenfe2ii2^. Fig.|4] summarizes the quali- 
tative agreement between temperature dependence of the 
dielectric constant (b) and the inverse of the phonon fre- 
quency (c). Remarkably, the measured changes of the 
dielectric constant correspond very well to a decrease of 
the phonon frequency as described by the LST equation. 
We note however that in order to have a quantitative 
agreement, one needs to add the contributions from the 
all the relevant symmetry allowed TO modes. 

In addition to this, the data imply a large anomalous 
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FIG. 4: Comparison between the behaviour of the phonon 
mode and the dielectric constant, (a) Temperature depen- 
dence of the phonon frequency responsible of the anomaly 
measured in the dielectric constant, (b) Dielectric constant 
measured perpendicular to the c axis, (c) Inverse of the 
phonon frequency as a function of temperature, (d) Line- 
width of the phonon mode associated to the dielectric con- 
stant anomalous behaviour. 



region close to the magnetic transition another decaying 
channel is available for the transverse optical mode. In 
the vicinity of the Neel temperature the phonon broad- 
ening deviates much stronger from the expected behavior 
suggesting that a phonon-magnon coupling channel be- 
comes active at the magnetic ordering temperature. This 
strongly suggests that the additional decaying channel re- 
ducing the phonon lifetime is not directly coupled to the 
magnetic order parameter but rather to its fluctuations. 
In conclusion, we have demonstrated magnetodielec- 
tric coupling in TbFca (603)4 using capacitance measure- 
ments and revealed the correlation between spin-phonon 
coupling and dielectric constant. We have shown that 
magnetodielectric coupling cannot be explained by the 
magnetostrictive effects but it occurs via shifting of the 
associated optical phonon modes. The magnetic field in- 
duced phase transition the magnetic field induced phase 
transition is associated to sharp changes of not only the 
magnetic but also dielectric and structural degrees of 
freedom, i.e. the spin-flop of Fe spins is accompanied 
not only by full polarization of the Tb moments but also 
to a significant distortion probably of the FeOe octaedra 
and a jump in the dielectric constant. 
Acknowledgments 



broadening of the phonon mode at the magnetic transi- 
tion itself. This is shown in Fig.|4l[d) where the phonon 
line-width as a function of temperature across the an- 
tiferromagnetic ordering temperature is presented. The 
black dashed line describes the expected decrease of the 
phonon lifetime as a function of temperature as if it was 
purely due to the usual phonon anharmonicity. The (red) 
dashed area in Fig.|4l[d) reveals that in the temperature 
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